EAULIQ: THE NEXT GENERATION 


by David A. Randall and Laura D. Fowler 



University 


DEPARTMENT OF 
ATMOSPHERIC SCIENCE 


PAPER NO. 



EAULIQ: THE NEXT GENERATION 


by 

David A. Randall and Laura D. Fowler 


Research supported by the 
National Science Foundation 
under Grant number ATM-9812384, 
by the U.S. Department of Energy 
under Grant number DE-FG03-95ER61968, 
and by the National Aeronautics and Space Administration 
under Grant Number NAG 1-1266. 


Department of Atmospheric Science 
Colorado State University 
Fort Collins, CO 


January 1999 


Atmospheric Science Paper No. 673 




Abstract 


This report summarizes the design of a new version of the stratiform cloud 
parameterization called Eauliq; the new version is called Eauliq NG. The key features of Eauliq 
NG are: 

• a prognostic fractional area covered by stratiform cloudiness, following the approach 
developed by M. Tiedtke for use in the ECMWF model; 

• separate prognostic thermodynamic variables for the clear and cloudy portions of each 
grid cell; 

• separate vertical velocities for the clear and cloudy portions of each grid cell, allowing 
the model to represent some aspects of observed mesoscale circulations; 

• cumulus entrainment from both the clear and cloudy portions of a grid cell, and 
cumulus detrainment into the cloudy portion only; 

• the effects of the cumulus-induced subsidence in the cloudy portion of a grid cell on the 
cloud water and ice there. 

In this paper we present the mathematical framework of Eauliq NG; a discussion of cumulus 
effects; a new parameterization of lateral mass exchanges between clear and cloudy regions; and a 
theory to determine the mesoscale mass circulation, based on the hypothesis that the stratiform 
clouds remain neutrally buoyant through time and that the mesoscale circulations are the 
mechanism which makes this possible. An appendix also discusses some time-differencing 
methods. 
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1 . Introduction 

What is a stratiform cloud? Virtually all stratiform clouds contain buoyancy-driven 
turbulence (e.g. Lilly, 1968). On the other hand, many convective clouds extend upward from and/ 
or grow upward into stratiform cloud layers. For purposes of this paper, we define stratiform 
clouds to be clouds which are, in an area-averaged sense, neutrally buoyant with respect to their 
environments, level-by-level throughout their vertical extents. In contrast, then, “convective” 
clouds are those which are either positively or negatively buoyant at some levels, in an area- 
averaged sense. For example, cumulus or cumulonimbus clouds are positively buoyant through a 
large fraction of their vertical extents, and they may also contain regions of negatively buoyant 
convective downdrafts. 

Many stratiform clouds are generated by convective clouds. Each day, tens of thousands of 
cumulonimbus clouds inject enormous quantities of boundary-layer air into the upper troposphere 
and lower stratosphere (e.g., Riehl and Malkus, 1958). The detrained air forms horizontally 
extensive and deep “anvil” clouds which contribute as much as 40% of the total precipitation that 
falls from the convective systems. Observations (e.g., Webster and Stephens, 1980) show that 
cirrus outflows from tropical convection can extend for many hundreds or even thousands of 
kilometers downstream from the convective disturbance that generates them. In such cases, the 
local time rate of change and advection terms of the conservation equation for the large-scale 
average cirrus ice water concentration must be comparable to the source and sink terms, so that a 
prognostic approach including advective effects is necessary for accurate predictions of the 
cloudiness. 

The anvil clouds contain mesoscale circulations (e.g., Houze, 1982) and also small-scale 
moist convective circulations, which influence the evolution of the convective systems and the 
large-scale circulations in which they develop. Houze (1982) showed that the total heating 
associated with a mature convective cloud system peaks in the upper troposphere and is quite 
small in the lower troposphere, whereas the convective heating alone is weaker in the upper 
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troposphere and stronger in the lower troposphere. Similar conclusions have been reached by 
Johnson and Young (1983), Johnson (1984), and Cheng and Yanai (1989). The heating and drying 
due to the mesoscale vertical motions are comparable in magnitude to the latent heating and 
precipitation drying associated with the anvil. The mesoscale heating and drying are not 
negligible, but neither are they as strong as the convective heating and drying. Rutledge (1986) 
and Rutledge and Houze (1987) have presented evidence that there are no fundamental 
differences between the stratiform anvils of midlatitude and tropical convective cloud systems. 
Houze (1982) has argued that mesoscale vertical motions associated with stratiform anvil clouds 
make important contributions to the large-scale heat and moisture budgets of convective cloud 
systems, and he has attempted to determine the typical magnitudes of these effects from 
observations. On the other hand, Cheng and Yanai (1989) have presented evidence that the effects 
of mesoscale condensation and evaporation are considerably larger than those of mesoscale 
vertical drafts. Certainly there is sufficient motivation to develop an accurate physically based 
parameterization of the effects of mesoscale vertical motions and mesoscale condensation and 
evaporation on the large-scale heat and moisture budgets. 

Anvil and cirrus cloud systems are generally too small to be explicitly represented in 
large-scale models; this leads us to the problem of “subgrid cloud amount,” which can be defined 
in general terms as the statistical distribution of cloud water and cloud ice on the subgrid scale. 
There are actually three reasons why we are interested in fractional cloudiness. The most obvious 
is that the small scale distributions of liquid water and ice can have a strong effect on the transfer 
of radiation. A second reason is that fractional cloudiness matters for cloud microphysics. The 
importance of this for realistic simulation of large-scale cloudiness has been emphasized by 
Fowler et al. (1996). As an example, the Colorado State University General Circulation Model 
(CSU GCM) described by Fowler et al. (1996) predicts the grid-cell averaged mixing ratios of 
liquid water and ice. Microphysical processes, such as conversion of cloud water to rain water, are 
local processes, and so they must be formulated in terms of the local concentrations of cloud 
water and cloud ice, rather than the predicted large-scale average concentrations. The local 
concentrations are essentially equal to the large-scale average concentrations divided by the cloud 
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amount. This problem has been discussed by Bechtold et al. (1993). A third reason for interest in 
fractional cloudiness is that it affects the dynamics of the convective clouds that produce much of 
the condensation. The effects of fractional cloudiness on cloud dynamics were discussed by 
Randall (1987). 

The earliest cloud amount parameterizations (e.g. Smagorinski, 1960) simply related the 
fractional cloudiness to the large-scale relative humidity i.e. 

f = f(RH). (1) 

Xu and Krueger (1991) tested this and other simple cloud amount parameterizations by using a 
Cloud System Model (CSM) 1 , and found that they do not work very well. A similar conclusion 
was reached by Xu and Randall (1996 a, b). 

In a very influential paper, Sundqvist (1978) proposed a stratiform cloud parameterization 
for large-scale models, in which the large-scale average cloud water mixing ratio was introduced 
as a prognostic variable, and simple parameterizations of microphysical processes were used to 
represent the sources and sinks of cloud water. The stratiform cloud parameterization recently 
incorporated into the Colorado State University General Circulation Model (CSU GCM), as 
described by Fowler et al. (1996), can be regarded as a recent attempt to follow the trail blazed by 
Sundqvist, using a somewhat more modem microphysics parameterization. Sundqvist related the 
cloud amount to the large-scale relative humidity, according to 


*• We use the term “Cloud System Model” to denote a model with sufficiendy high spatial resolution to 
resolve individual cloud elements, and integrated with a domain size large enough to encompass many indi- 
vidual clouds, and with a time domain long enough to include many cloud life cycles. A second term some- 
times used is “Cumulus Ensemble Model,” but this seems inappropriate because the models can be applied 
to stratiform clouds. A third term sometimes used to refer to these models is “Cloud Resolving Models.” 
This name is not very satisfactory because it makes no reference to the large domain size and long integra- 
tion time characteristic of such models. 
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/ = Max 


( 2 ) 


f RH-(RH) 0 J 

I 1-Wo* J’ 

where (RH) 0 is a “threshold” relative humidity below which the cloud amount is assumed to be 
zero. Sundqvist further assumed that when stratiform clouds exist, the relative humidity in the 
cloud-free portion of the grid cell remains constant at (RH ) 0 , while the relative humidity in the 

saturated portion of the grid cell is equal to the saturation value there. If the area-averaged relative 
humidity is predicted, then this assumption permits diagnostic determination of the cloud fraction. 
As explained later, in Eauliq NG we do not assume that the relative humidity in the clear portion 
of the grid cell is a constant. 

A key assumption of Sundqvist (1978), which is not explicitly stated in the paper, is that 
the potential temperature is horizontally uniform throughout the grid cell. This is roughly 
consistent with the definition of stratiform clouds given above, i.e. stratiform clouds are neutrally 
buoyant (in an area-averaged sense) with respect to their environments. The parameterization 
described in the present paper makes use of this “neutral buoyancy” assumption, but a key point is 
that we specify a physical mechanism which maintains this neutral buoyancy over time. 

A different approach was suggested by Sommeria and Deardorff (1977), who proposed a 
subgrid cloudiness parameterization intended for use in high-resolution, cloud-resolving models. 
They assumed that a pair of moist conservative variables, such as liquid water potential 
temperature, 0/, and total mixing ratio, q t , undergo subgrid-scale fluctuations and have a joint 
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Gaussian probability density function (pdf), as sketched in Fig. 1 . The shaded region of the sketch 



Figure 1 : Sketch illustrating the Sommeria-Bearterff approach to cloud amount 

parameterization. The shaded portion of the figure is unsaturated, and the white portion 
is saturated. The @va! represents an isollne of the joint probability distribution of 6/ and 

<7f 

represents subsaturated air, and the white region represents saturated air. By integrating the joint 
pdf over the saturated region, Sommeria and Deardorff were able to determine the cloud amount. 
In order for this approach to be workable, it is necessary to know the joint pdf of the two moist 
conservative variables. Le Treut and Li (1988), Smith (1990), and Ricard and Royer (1993) have 
followed simplified versions of this approach in large-scale models. Le Treut and Li (1988) and 
Smith (1990) used assumed (rather than predicted) pdfs for the moist conservative variables. In 
principle, these pdfs should be determined by small-scale and/or mesoscale dynamical processes, 
such as convective turbulence within the stratiform cloud (Randall 1987; Randall et al. 1992). 
Ricard and Royer (1993) implemented their parameterization in the context of “level 2” subgrid 
scale turbulence parameterization (Yamada and Mellor 1979), using 20 layers to represent the 
vertical structure of the atmosphere. Xu and Randall (1996 a) used a CSM to evaluate the 
applicability of pdf-based parameterizations to the simulation of large-scale circulations. They 
found that the coefficients of such parameterizations are cloud-regime-dependent Nevertheless, 
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pdf-based parameterizations certainly have the potential for further improvement, and in fact the 
parameterization presented in later this paper can be interpreted as a pdf-based scheme. 

Albrecht (1981) considered the gradual evaporation of cloud water produced by 
detrainment from shallow cumulus clouds. He proposed a parameterization of the form 


^^cld _ ^ v ^clr 


(3) 


where (q,) ,, is the total in-cloud mixing ratio (vapor plus liquid), q * is the large-scale saturation 

mixing ratio, and (# v ) clr is the vapor mixing ratio of the clear air. Albrecht assumed that the 

cloudy and clear portions of the grid cell have the same temperature. His parameterization was 
based on a simple model in which the detrained cloudy air was assumed to “relax” back towards 
the mean-state mixing ratio. A simple convective cloud model was used to estimate (^) cld • The 

parameterization invoked microphysical processes such as the evaporation of cloud water, but did 
not explicitly parameterize the microphysics. S. A. Klein (personal communication, 1997) has 
pointed out that (3) is equivalent to 

f = (4) 

+ ')«• 

This means that the cloud amount increases as the in-cloud liquid water mixing ratio increases, 
and as the large-scale relative humidity increases. Note, however, that (4) can give / -» 1 even 

for RH < 1 . Xu and Randall (1996 b) proposed a semi-empirical cloud parameterization which is 
somewhat similar in spirit to (4): 
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(5) 


/ 


RH P \ 


1 - 


exp[-g(g,) cld J 

[(1 -RH)TJ 


They used a CSM to evaluate the parameters used in (5), for two different cloud regimes. 

Tiedtke (1993) developed a cloud amount parameterization for use in the ECMWF model. 
He introduced a prognostic equation for the cloud amount in addition to prognostic equations for 
the mass of cloud water. This prognostic equation was a major advance of Tiedtke’s approach; it 
simply expresses the conservation of mass for the cloudy air, i.e. it is essentially a continuity 
equation. One of the strengths of Tiedtke’s parameterization is that convection acts as a source of 
stratiform cloud water and cloud amount; this idea was inherent already in the cumulus 
parameterization of Arakawa and Schubert (1974). Tiedtke also included the effects of the 
advection of condensed water and cloud amount by the large-scale circulation, and he 
parameterized the effects of microphysical and turbulent processes on the clouds and the large- 
scale thermodynamic state. In accord with the definition of stratiform cloudiness given at the 
beginning of this paper, Tiedtke assumed (tacitly) that the temperature is horizontally uniform 
across the clear and cloudy sub-regions of each grid cell, but he did not explain how this uniform 
temperature is maintained. 

Fowler et al. (1996) and Fowler and Randall (1996 a, b) developed and tested a bulk cloud 
microphysics parameterization called “Eauliq,” based on the work of Rutledge and Hobbs (1983) 
and Lin et al. (1983). Eauliq included representations of the microphysical processes responsible 
for the formation and dissipation of both water and ice clouds. Eauliq included five prognostic 
variables representing the mixing ratios of water vapor, cloud water, cloud ice, rain, and snow. 
Graupel and hail were neglected. Cloud water and cloud ice were permitted to form through large- 
scale condensation and deposition processes. Rain and snow were assumed to be produced 
through autoconversion of cloud water and cloud ice. Rain drops falling through clouds were 
assumed to grow by collecting cloud water, and falling snow was assumed to collect both cloud 
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water and cloud ice. These collection processes were formulated using the continuous collection 
equation. Evaporation of cloud water, cloud ice, rain, and snow were allowed in subsaturated 
layers. Melting and freezing were considered. Fowler et al. also included a coupling between 
convective clouds and stratiform anvils through the detrainment of cloud water and cloud ice at 
the tops of cumulus towers. Interactive cloud optical properties provide the link between the cloud 
microphysics and radiation parameterizations; the optical depths and infrared emissivities of 
large-scale stratiform clouds were parameterized in terms of the cloud water and cloud ice paths. 
Perhaps the most serious weakness of Eauliq is that it does not include a parameterization of 
cloud amount; Fowler et al. (1996) simply assumed that the cloud amount was either zero (i.e. no 
cloud in a grid cell) or one (i.e. uniform cloud throughout a grid cell). 

The purpose of this paper is to outline a generalization of Eauliq, called “Eauliq NG,” 
which includes a variable cloud amount. We determine the cloud amount prognostically using an 
equation similar to that proposed by Tiedtke (1993). We also diagnose the thermodynamic 
properties of the clear and cloudy portions of the grid cell, using the prognostic cloud amount, the 
prognostic mean-state thermodynamic variables, and the prognostic differences in the 
thermodynamic variables between the clear and cloudy portions of the grid cell. We separately 
determine the vertical motions in the clear and cloudy portions of the grid cell, by requiring that 
the differences in vertical motion act to maintain neutral buoyancy of the stratiform clouds. These 
same vertical motion differences have additional consequences, of course, which are taken into 
account in our parameterization. 

Section 2 of this paper lays out the conceptual and mathematical framework of Eauliq NG. 
Section 3 discusses the role of cumulus processes. Section 4 outlines the algorithm used to 
determine the separate thermodynamic properties of the clear and cloudy sub-regions of each grid 
cell. Section 5 describes our method to determine the separate vertical motions in the clear and 
cloudy sub-regions. Section 6 discusses the parameterized lateral mass exchanges between the 
clear and cloudy sub-regions. Section 7 describes the parameterized microphysical processes. 
Section 8 gives a summary and conclusions. 
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2. Framework 

2. 1 Sub-regions 


We divide the horizontal 2 domain of a grid cell into three sub-regions, as shown in the 
sketch below. These are the clear sub-region, denoted by “clr,” the sub-region filled with 
stratiform cloud, denoted by “cld,” and the cumulus sub-region, denoted by “cu.” The areas 

occupied by the three sub-regions are denoted by Al clr, A l c id, and A l C u, respectively. The total 
area of the grid cell is 


A — A clr + A cld + A cu • 


( 6 ) 


Here the superscript i denotes the grid cell under consideration, and subscripts are used to denote 


a sub-region. We assume that Al is independent of both time and height. 



Figure 2: Schematic showing the horizontal cross section of a grid cell containing a clear 

region (clr), a stratiform-cloudy region (cld), and a cumulus region (cu). 


2 We assume that the model has a vertical resolution high enough to capture at least crudely the vertical dis- 
tribution of the cloudiness; we therefore ignore the possibility of vertically subgrid-scale clouds. 
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By definition, the cld sub-region contains stratiform cloud water and/or cloud ice, and by 
definition the clear region contains neither. In other words, cloud water and cloud ice are assumed 

to occupy the same fractional area, i.e. A 1 c u • We allow the possibility that rain and snow can exist 
in either cloudy or cloud-free portions of a grid cell. Obviously rain and snow must originate in 
cloudy regions, but they can fall into clear regions. The fractional areas occupied by rain and snow 
are discussed later. 

In Fig. 2, the cumulus sub-region is sketched as if it were comparable in size to the others, 
but this is merely for convenience in making the drawing; in reality we expect the fractional area 
occupied by the cumulus clouds to be very small in all cases. The “cu” sub-region represents an 
ensemble of cumulus clouds, which may be further broken down into subensembles, following 
the “spectral” approach of Arakawa and Schubert (1974). Here we forego the additional 
notational complexity that would be required to explicitly represent the subensembles, but we do 
in fact follow the spectral approach. 

Following the approach of Margolin et al. (1997), we allow each of the three sub-regions 
(“clr,” “cld,” and “cu”) to exchange mass laterally with the other two. This mass exchange can 
occur “inside” the grid cell under consideration, and in addition, each sub-region in grid cell i can 
exchange mass with neighboring grid cells. We denote lateral mass exchanges inside the cell 
(hereafter “intra-cell”) by E, and lateral mass exchanges with neighboring cells (hereafter “inter- 
cell”) by F. 

Consider an arbitrary intensive variable h , and let S h denote the source or sink of h , 
which can include the effects of small-scale turbulence and radiation, as well as microphysical 
processes 3 . We can write the following budget equations for grid box i : 


3- For now we use h as a generic intensive scalar; later we use the same symbol to denote the generalized 
moist static energy, which is of course an intensive scalar. 
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( 7 ) 


bt c\iA l ck) = E l c\d,ckh l cld + E* cu,ckh cu - ( E clr, cld + E clr,cu)^ ch- 
at 

- Y^E 1 ’ 1 clr h 1 ’ 1 ck-^mwckh clr^clr) + (^) i c lr^ IcLr ’ 


-^-(AW^cld-^cld) = ^clr.cld^clr + F^cld^ cu - (E cld, clr + E cld,cu)^ cld 
Ot 


-2/' cld ^’ 1 cld - S-(m‘ w l cldh cld^cld) + (*S'/,) i c ld v4 * cld ’ 


dz 


( 8 ) 


h' cuA' cu) = E i clr.cu^clr + F^cld.cu^ cld - (E cu,clr + E cu,cld)^ cu 
ot 

- jV ’ r cu h u f cu - ^(mJcuhcuA^u) + (Vcu^'cu . 


(9) 


Here m is the mass of dry air, per unit area, in grid cell i . The E ’s and F ’s have dimensions of 
mass per unit time. The terms involving the E ’s and F’s are discussed in the next two sections. 
When we add (7), (8), and (9), all of the E terms cancel out, but the F terms survive. We obtain 

st mhA) - 

( 10 ) 

+ A\S h ) 1 , 




i, i ft,t 

clr/* 


clr + F 1 ’ 1 cld h U 1 cld + F 1 ’ 1 cu h '’ * cu) - -^JiA'mwh 


where 


A l h l = h l c ]iA l c]i + h l c ldA l c\d + h' cu A cu, 


( 11 ) 


A* mwh = m\v ckh clr^clr + m‘w cldh cldA cld + m' w cuh cu-<4 cu, 


( 12 ) 


and 
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As h )‘ = (Sj) , cll ^'cl I + (S A )' cld ^'dd+(S 4 )' cu ^'cu. 


( 13 ) 


The continuity equations corresponding to (7)-(10) can be obtained by setting h= 1 and 
S h = 0: 


— (tw A clr) — E cld.clr + E cu.clr ~ (E clr.cld E clr.cu) ~ clr — w clr-^ clr), (14) 


dz 


— (m A cld) — E clr.cld + E cu.cld~ (E cld.clr + E cld.cu) ~ cld — w cld-^ cld) , (15) 


— (m A cu) — E clr.cu + E cld.cu ~ (E cu,clr + E cu.cld) ~ "^.E M, cu ^ cu)) (16) 


J^( m l A l ) = ^(^’‘ck + F 1 ’ 1 cld + E 1 ’ 1 cu )--g-(m l w A 1 ). 


(17) 


In (17), 


i i j i f i i 1 

Aw — w cliA clr w cld-^ cld w cu-4 cu » 


(18) 


where w is the area-averaged vertical velocity. Eqs. (14)-(16) govern the time change of the mass 
or area within each sub-region. Eq. (17) is the continuity equation for the whole grid cell. By 
combining (7)-(9) with (14)-(16), we can derive “advective forms” of the three budget equations: 
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( 19 ) 


m 


‘ A* ckS-ih* ck) ~ E‘cld,c\i(h cld — h clr) + -£ cu,clr(^ cu~ h clr) 
ot 


-X F ' ,r clr(^’ I clr-/l'clr) 

r 



A‘ c\d^~(h‘ c\d) ~ E‘ck,c\d(h clr — h cld) + E cu,cld(A cu _ h cld) 
Ot 


^ F *' 1 cld(^’ * cld - h 1 cld) 


— m A cldw c\&Q^h cld) + (S^) cld 


( 20 ) 


m A cir~r(^ cu) — E> clr,cu(^ clr — h cu) + E cld,cu(^ cld - h cu) 
ot 

-JV’ 1 ’ cu(h 1 ’ 1 cu-h'cu) 


( 21 ) 


m A cuW ca ~Q^J l cu ) + (^h) CU A cu 


2.2 Lateral exchanges of mass between subregions 

First, we discuss the mass exchanges between cumulus clouds and the clear and 

stratiform-cloudy portions of the grid cell, f'cu.clr represents the flow of mass from cumulus 
clouds into the clear portion of the grid cell; here we adopt the convention that the first subscript 
denotes the sub-region of origin, and the second denotes the destination sub-region. Similarly, 

£'c U) cld represents the source of stratiform cloudy air due to detrainment of air from cumuli into 
the cloudy part of the grid cell. We assume that air detrained by cumuli is always cloudy, and so 
always enters the stratiform-cloud sub-region of the grid cell. This means that 


E cu,clr — 0 • 


( 22 ) 
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For notational clarity, we write 


Z>Wcu,cld>0. (23) 

£'clr,cu represents the entrainment of air from the clear part of the box into cumulus clouds. 

Similarly, £ ,, c ld,cu represents the entrainment of air from the cloudy part of the box into the 
cumulus clouds. We note that 


E clr,cu — 0 , and E cld,cu — 0 . (24) 

Now consider the mass exchanges between the clear and stratiform cloudy portions of the 

grid cell. i^cld.clr represents air “moving” from the cloudy area to the clear area. Of course 

£*cld,clr amounts to evaporation of the cloud, which is not really motion in the usual sense at all. 

Similarly, £ i c ir,cld represents the transformation of clear air into stratiform cloud air. Again, this 
is not really “motion;” instead it represents the effects of processes that increase the relative 
humidity in the clear part of the box, or some portion of it, so as to convert the clear air into 
cloudy air. Note that, with these definitions, 


£*clr,cld > 0 , and £*dd,clr ^ 0 . (25) 

In particular, £*clr,cld is not equal to minus £ I cld ) clr (unless they both happen to be zero). The two 
mass exchange processes can occur independently and simultaneously within each grid cell. 

When we apply (19) to the cloud water mixing ratio, q c , which is of course equal to zero 
in the clear portion of the grid cell, we obtain a very simple result: 
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Note that here q 1 denotes the value of q . in the cloudy portion of the grid cell; we omit the 

subscript “cld” on the grounds that it would be redundant. According to (26), the flux of cloud 

water from the cloudy to the clear portion of the cell, denoted by £*cld,clr(tf c ‘) , is balanced by a 

sink of cloud water (presumably due to evaporation), which prevents any cloud water from 
accumulating in the clear portion of the grid cell. We can use (26) to diagnose the rate of 
evaporation in the clear-portion of the grid cell; not only is this needed for such diagnostic 
purposes, it also represents a source of water vapor in the clear portion of the cell, and it is 
associated with evaporative cooling there. Similarly, the cloud ice mixing ratio satisfies 

0 = E cld,clr(^ , j ) (6'_ ) . ^ clr • (27) 

1 * i clr 

Intra-cell fluxes between sub-regions of the same type (“sub-sub-regions”) are neglected in 
this paper, so that, for instance, flows between two clear portions of the same grid cell are 
assumed to have no effect on quantities of interest. In effect we assume that all subregions of the 
same type, within a given cell, have identical properties, so that exchanges among them are 
irrelevant. 

2.3 Inter-cell exchanges of mass 

We assume that when air flows across cell walls, between neighboring grid cells, it always 
moves between sub-regions of like type, so that, for example, air can travel from the cloudy sub- 
regions of grid box i to the cloudy sub-regions of neighboring grid box f , but not from the cloudy 
sub-region of grid box i to the clear sub-region of neighboring grid box /' . The rationale is that it 
would be quite unlikely for the boundary of a sub-region, e.g. a stratiform cloud, to coincide 
exactly with the wall of a grid cell. 



To minimize the number of symbols, we adopt the notation F h 1 c lr to denote the flow of 
mass outward from grid box i to neighboring grid box t , in this case between the clear sub- 

regions of each. Corresponding conventions are used with F c ld and F ’ C u • This is why the F 
terms appear with minus signs on the right-hand sides of (7) - (9). The F s can have either sign. 

1^' clr represents the net flow of air from the clear region in one cell to the clear 
f 

j ji 

regions of the surrounding cells. Similarly, ][]F ’ c ld represents the flow of air from the 

r 

stratiform cloud region in one cell to the stratiform cloud regions of the surrounding cells, and 

z z 1 

^T\F ’ cu represents the flow of air from the cumulus region in one cell to the cumulus regions of 

i' 

the surrounding cells. Summations such as ^F *’ 1 c lr h ' clr represent exchanges between grid 

f 

* i t 

box i and all neighboring grid boxes. The “hat” symbol, as in h ’ , denotes an interpolated value 
on a cell wall. 

For simplicity, we currently neglect any variations of the horizontal velocity between the 
cloudy and clear regions of the grid cells, although we realize that this must be addressed in the 
future. We assume that 


, i, f i, e .i, i 

fn v nl clr, 


.* j 1 f f 

where m ’ is an interpolated mass variable, defined on the cell wall, v ’ n is the outward normal 

velocity component along the cell wall, and / ’ clr is the distance occupied by clear air, along the 
cell wall. Similarly, we assume that 
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jJ, i . i, t i, f ji, f , , 

F cld = m v nl eld, 


( 29 ) 


and 



„ i, ‘ i, i ,h i 

m v n l cu ■ 


(30) 


We require that 


/ clr + l cld + * cu — * 5 


where /*’* is the total length of the cell wall. We assume that 


. j, i 

ji, v zM clr 


/’clr = / 




-i, i 

j, r _ _ ,1, jM cld 


7*J * 

/ cld = / 


/, z 


„ V * 1 

f z, f _ ji , / cu 
* cu — l 7“T" ? 


(31) 


(32) 


(33) 


(34) 


where 1 clr, 1 cld, and ^4*’ 1 cu are suitably interpolated “edge” values of the corresponding 
areas. In view of (31), we must require that 


-i, f «/, v , - i * i, i 
A clr + A cld + A cu — A 


(35) 


This requirement is met by either upstream or centered interpolation of the areas, and also by any 
linear combination of upstream and centered interpolations. Substitution of (32) - (34) into (28) - 
(30) gives 
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( 36 ) 


y, ‘ 

rJ, i * i, f i, f ,i, FA clr 
F clr = m v n l — 

A . * 5 * 


- J, l 

-,i,V _ . F i, F ji,fA cld 


Ft, t A *5 * »> * /*! 

cld = m v n / 




( 37 ) 


• , A 1 ’ 1 

t CU 

r cu tw v «/ — t-t- 

y» i 


( 38 ) 


Substitution of (36)-(38) into (17) leads to 


d , / y> *\ x -1 

rt mA >-? 


m v nl 


{ ~i,f *i,F ~i, F V 
.4 clr + .4 cld + -4 cu 




yj 


~g£ mwA > 


i, T i, T j, f 5 / i—F i,f 
m v n l --^(mwA ). 


( 39 ) 


To obtain the second equality in (39), we have used (35). Eq. (39) is just what we would expect 
from the large-scale point of view. On the other hand, Eqs. (36)-(38) imply a kind of “diffusion” 
of the mean state properties across cell walls. To see this, substitute (36)-(38) into (10), to obtain 


£ (F 1 ’ 'dr ti' ^clr + F U f dd ti' 'cld + F U f CU ti' 'cu) 


( 40 ) 


-2 


f 

^ „ i, F i, F ,i, CA ’ clr fU? . , * i, F F ,i, F A cld £*> f 

m v " 4 ■■ 


, , A 


v 


yi, r A clr?M . , A i, f i,f ji,fA cld fh* % , , a i, T i, f , 1 , M cu f U i 
nl — —h clr + m v n l — rrrh cld + m v n l — —h cu 

\ l > 1 *A 1 ’ 1 £ 1 ’ 1 




^4 clr + A cld + ^4 cu 

\U f 


h 1 ’ 1 -£(?'' (h) 
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where 


\i, 1 . , 

3( ? 


-I,/ . - 

Z A z + f i + f 7 z + zM cld ?h 1 , , 
m V n l cld 

A *' 1 

r A 


*J, z 


Z A /, f z, f 7 i,fA’ cuth? 3. z t} A K , a*7q~V 
m v «/ — —A cu-^(m w/zi) + ^(^) 

z’ 1 .4 ’ 


^ ’ clr + 4 cld + v4 


S * h f z, f 7 z, z 1 

m v nl 


S ’ 1 ' ^ 


cu 




?Z, z 

A 


( 41 ) 




Z 


5> 4 V’ V’ V f + J^G?' l \h)~^-(mwhA l ) + A\S h , ) , 

f l 


Z, Z 


s-yi, ? / 1 \ *h? hi jhi4 clr ,fhi ?z, z 

G (h) = -m v /j/ — — (/* clr — A ) 

^ z, z 
4 


A z, z " 1 z, z * 1 7 z, zM ’ cld^ f z, r *i, i\ 

-m V n l rrr(h cld -h ) 

* z, z 


(42) 


- i, ? i,f ,i,?A’ cu ,£*»*’ 

-m v «/ — — (« cu h ) 

- Z, z 

4 


represents a diffusive flux that arises from the indicated differences in h . To obtain the third 
equality in (40) we have again used (35). 

2.4 Vertical mass fluxes 

It is useful to express all vertical mass exchanges in terms of three mass fluxes: 
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( 43 ) 


i i i 1 i i i i 

the large-scale mass flux, A M = A c l iM c lr + A c id A/ c ld + A C u M C u ; 


the cumulus mass flux, M * c = — — (A/ C u - M ) ; and 


(44) 


.i J 


the mesoscale mass flux, mcso —(M 1 cld - M* c\i) ■ 


A l A l 


(45) 


Here 


A/ 7 clr = rnw clr; 


(46) 


M 1 cld = m l w l dd ; 


(47) 


M 1 cu = mw l ea- 


rn 


Straightforward algebraic manipulations, which are summarized in Appendix A, lead to the 
following results: 


A tA ^‘ddA i ch /1 A i A . \ A c kA C u, t A 


A'clrM clr = .4'clrM - ° ld - (H/dd - M* clr) — -(A^cu - A/ C lr) , (49) 


1 A t~J , A c\dA c ]r A *A . \ cld A cu , , j 


A 1 c \dM cld = A'cUM + Cld - (A^dd ~ M* clr) ^^^r-^Wcu - A/dd) , (50) 


mwh‘ = M 1 h l + M l meso(h t cld-h‘c]i) + M I c(h l cu-h). 


(51) 
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2.5 Summary of framework 

Before making further simplifying assumptions and approximations, we summarize the 
basic equations of our model, making use of the various results obtained in the preceding sub- 
sections. The continuity equations can be written as: 


-^-(m^clr) — ■E , cld,clr + E cu,clr — (E clr, cld + E clr.cu) - 
ot 

A f 

z, f A,? A clr d , A i %jf l \ 

m v 9 nl —-^4 clrMclr) , 

V’ * oz 
f A 


(52) 


^(jn A* cld) “ E clr, old + E cu,cld — (E cld,clr + E cld,cu) ~ 
Ot 


(53) 




z, r z, r 7 i, i A cid d f A i _ %Jf l , ^ 
m v n l — — (^ c\&M c ld) , 

-J, 1 OZ 

r A 


^-(m A' cu) “ E clr,cu + E cld,cu ~ (E l CU, clr + E cu,cld) — 
ot 


(54) 


* ; r A cu d 


Z * z, f z, z 7 z, z A CU O ( J V 

m v n l — - —(4 cuM cu , 

dz 

f A 


d / i A - U f U f iU f d , z_* A 

— (mA) = ~2_ J m v n l (w w A). 


(55) 


The conservations equations for h can be written as: 


— (m^clr-^clr) — E cld, clr/* cld + E cu,clr/* cu — ( E clr, cld + E clr,cu)/* clr 
ot 




- £ m‘- 'V' *'• 'dr - ^'drJu'drA'cl.) + (S 4 )' cl /c . 

f A ’ 


(56) 
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— (/w h cld^4 cld) ~ E clr,cld h clr + E cu,cld h cu — ( E cld,clr + E cld,cu)^ cld 
ot 


, I 


- JV’ V’ V’ A 1 ’ ' dd - ^f^ddMdd/i'cld] + (^) , cld ^ / cld , 

f A ’ 


( 57 ) 


iitnh' cuA cu) — £ clr,cu^ clr + £ cld,cu^ cld — (£ cu,clr + £ cu,cld)^ cu 
Ot 


~.i, i 


Z «i, f f ,i,i'A’ cut d , ./ >/ ,i x , / C V 

m v nl — :rp-/j cu-^cuMcuAcu) + (^) cu ^cu 5 

/’ ^4 ’ 


( 58 ) 


d , iji , is y 1 * i, /’ i, J 1 »/, f tU ? , f / v \ 
— (m h A) = - > m v «/ h + G (h) 
ot , 


— { ^4 [A/* A + AA* meso(h cld — h clr) M* c(h cu — h )] } + A (5^) 


( 59 ) 


Our strategy for the use of these equations is explained later. 

3. The cumulus terms 

We assume that 


r i _ A clr r i 
& clr,cu - — —-£• » 

A 


( 60 ) 


„i _ A cld r,i 
E cld,cu “ —t, . 

A 


( 61 ) 


Here E l is the total rate of entrainment into the cumulus clouds, i.e. 


E — E clr,cu ■*" E cld,cu • 


( 62 ) 
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Eqs. (60) and (61) mean that the air entrained into cumuli is derived from the clear and cloudy 
subregions in proportion to the fractional areas covered by those subregions. This assumption is 
made for simplicity and should be refined in the future. 

We also assume that the cumulus clouds detrain only into the cloudy portion of the grid 
box, on the grounds that the air leaving the cumuli is, after all, cloudy: 


E cu,clr — 0 . 


(63) 


To make our notation closer to that used in the earlier literature, we write 


E cu,cld = D . 


(64) 


Next, we assume, following Arakawa and Schubert (1974), that the fractional area covered 
by cumulus convection is very small compared with one, i.e. 


. ; i , 

A cu « A and A cu « A , 


(65) 


and also 


A*ck + A 1 d d = A i and A*’ '’clr + A*' '"eld 5 A *’ ' 


( 66 ) 


This assumption allows us to simplify our equations as follows: 


A h = h clr^4 clr + h cld-^ cld 5 


(67) 


on, 

A 1 


( 68 ) 
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F^cusO, 


(69) 


„U . i, f i, f ,i, t A l " , ,A r , f 

G (h) = -m v n l A ck(h ck-h ) 
-m v n l A cld(/i cld - « ) , 


(70) 


j * j i i i i i 

A ckM clr = -4 clr-A/ — A M meso A ckM , 


(71) 


i J i — ^ i i i 

A c\&M cld = A cld M + A M meso — A c\&M r 


(72) 


In addition, we can neglect the time-rate-of-change terms in (58) and (54). 

Using the various assumptions and definitions discussed above, the continuity equations 
(52)-(55) can be written as 


A‘ ck) — — v’ nl’ A ck + F cld,clr — F clr, cld 

ot “ 

i 

- |V dr M i ’ - A l M meso) + ^UmJ A* ck) - , 


(73) 


-^-(/n ^4 cld) — — ^ JW ’ v ’ nl A cld + E clr.cld — F cld.clr 
Ot ^ 


- h/cldM 1 + A* M meso) + Id) - F^ + D' , 


(74) 


0 = 


(75) 


d, * A v 1 * *» I > 1 ’ /M’ d sTZ* aK 
— (m v4 ) = 2," 1 v « Z A ) 


dz v 


(76) 
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The vertical mass flux divergence terms of (73)-(74) have been expanded to separately exhibit the 
contributions associated with the mean flow, with the mesoscale circulation, and with cumulus 
convection. Similarly, Eqs. (56) - (59) can be rewritten as: 


tt-(/w h clr -<4 clr) — E clcfclr^ cld — E clr,cld^ clr v n ^ 


1 


dt 


. i, f i, i 1 ,i, f A clrf ! > ? 

•• ft clr 


' J, i 

A 


d_ 

dz 


J^-^clr-W - A‘M meso)/2 + (S^) ^A clr 


(77) 


0/',. i .i J \ A clr u i , / 

+ dz[ Mc A chh clr J — r E h ck ’ 

^ A 


d , i , i A N i i r~<i i i * U ? U ? iU ? ^ cld L z ’ 1 

— (m h ckM dd) = E clr, cld h clr -E cld, clr h cld-^w v n l — 

dt f A ’ 

-J^(^ l cldA/ + meso)/l c idJ + (‘S , / J ) , c ld^ Ic,d 


cld 


+ MjA l Mh dd ) - ^-fE l h‘ dd + Ah' cu , 

7 ,4 


a_ 


(78) 


f J 


0 = E 


A clr , i , A cld , i 
— -h clr + -h dd 

V A 1 A 1 j 


-D l h l 


cu 


-A(A/ c , ^ , A < cu) + (5 a O , 'cu 


(79) 


t 

■ ^-{ A^A^h + M 1 meso (h cld - h clr)] } + [(^) c id^ cld + (*$/,) q^A clr] 
oz 


(80) 


--^[y4 M 1 c(h cu - h )] + (S^) CU A cu • 
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In (80), we have expanded the source/sink terms to show the clear, cloudy, and cumulus 
contributions separately. In each of (77), (78), and (80), we have written the cumulus terms at the 
end. Eqs. (75) and (79) are straightforward generalizations of corresponding equations proposed 

by Arakawa and Schubert (1974). 4 

As discussed by Arakawa and Schubert (1974), we can use the cumulus cloud model 
represented by (75) and (79) to write the cumulus terms of (80) in an alternative form, as follows: 


-(A* mV) ~ ~[A* Me (V cu - h 1 )] + ^cu^cu 


dt 


P » i j i i c)h ' i i 

= -^-(A l M c h l cu ) + A l M c ^ +Vh 7r / +A l cu(S l h ) cu 


D l h l cu~E l 


dz 

\ 


^ A clr , i A cld , i 
— -h clr + ~h cld 

{A* A 1 ; 


1 A 


- ( s h ) cu A cu 


+ A*Mc ~ +h i (E i -D i ) + A i cu(S i h ) cu 


= A 1 Me $■ + D^Vcu - V) + E l 
oz 


A clr,r» ,i s, . A cld ,ri ,i x 
~(h - h clr) + ~(h ~ h cld) 

.A 1 A 


(81) 


To obtain the second equality in (81), we have used (75) and (79). The entrainment terms on the 
last line of (81) did not appear in the corresponding equation of Arakawa and Schubert (1974), 
and represent a minor generalization of their results. These terms appear because of the 
differences between the cloudy and clear subregions. 

The cumulus terms of (74) can also be rewritten in a more useful form, by using (75): 


4 Note, however, that the definitions of E and D used by Arakawa and Schubert (1974) included a normal- 
ization by the grid-cell area; no such normalization is used here, which is why the grid-cell area appears in 
the third term of (75). 
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5 / / ,i \ d f i £ i * i \ cld , r\^ 

^-(m A cld) ~ ~-(M A dd) - E — — + D 
dt dz c A l 


i dM c 

= A + M cg z 


dz 
id A 1 dd 


D* + £ 

OZ 


/ x jr i a^\ ^ cl 

(M A ) — - 

J A 1 


A cld + D i 


idA cld 

~ M ‘Tz + 


f J 


1 - 


A dd 

A* J 


ID 1 . 


(82) 


This form agrees with that proposed by Tiedtke (1993) and derived heuristically by Randall 

(1995). It is useful because it shows clearly that the detrainment term does not try to drive .4* cld to 
values larger than one. Similarly, the cumulus terms of (73), (77) and (78) can be rewritten as 
follows: 


d , i J x d , , , i ,i s pi A clr 

of** A clr) ~ fa( M c A dir) ~E — 

A 


A dM c .u as idA ‘ cb Tn 1 + d ( u '■ A"K clr 
- A ‘ +m ‘Tz -\_ D+ d£ M c A \ J-T 


(83) 


idA clr jjA clr 

= M ^ — —V 

c dz A 


d (mV'clr^dr) ~ ^(M^'clr^'clr) - £ V c lr^ 
az A i 


dt 


,i dM c 


ti c\rA.‘ c\i-q^ * + M^ticUch) - [V + ^mJa^ 


(84) 


= ck^ dr) -D l h l c]T ^r , 


c dz 


i , i ^4 dr 
i 
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(85) 


^(mVcld^cld) ~ cUM'cld) - E‘h l c \d— + D‘h‘ on 

Ot OZ c 


= Add^dd^ + Mj^cld^dd) - |V + jL(M//)]A / dd^ + D l h 


i 8 ,, i 


i , d 


i A dd , ^,i , i 


cu 


— M c (A cld^ eld) D | h cu — h dd 


i , i A cld^ 


A 1 ) 


^ ^ . . 

It may seem surprising that cumulus detrainment influences —(m A c lr) and — (/w A c lrA c lr) , as 

indicated in (83) and (84), but the explanation is simple: detrainment decreases the clear area 
because it increases the area occupied by the stratiform cloud. In (83), the detrainment term 

decreases A 1 dr, but because this term is proportional to A l c \r it will never make A l c b < 0 . 

Similarly, the detrainment term of (82) tends to increase ^cid, but it will never make y^cld > 1 

because it is proportional to 1 - — c - . 


4. Parameterization of intra-cell mass exchanges 

4. 1 The relationship between cloud perimeter and cloud area 

The boundary between the cloudy air and the clear air is an obvious site for evaporation 
and/or sublimation to occur (hereafter we use the term “evaporation” as a shorthand for both 
evaporation and sublimation). The lateral mass exchanges between the clear and cloudy sub- 
regions of a cell occur across cloud lateral boundaries, so the rate of such mass exchange would 

be expected to vary in proportion to the cloud perimeter. From this point of view, E l c id, clr and 

A^clr.cld should depend in part on the amount of lateral “cloud edge,” or the length of cloud 
perimeter, that exists within the cell. 5 


29 



There are many possible relationships between the total cloud area and the total cloud 
perimeter. To illustrate this, we consider two possible ways in which the cloud area could 

2 

increase. As a first example, suppose that a stratiform cloud is divided into n identical square 

sub-clouds, each of area a = N 2 /n and each separated from the others by clear spaces of the 
same size as the clouds. This “checkerboard cloud” is illustrated in the left-hand panel of Fig. 3. 




Figure 3: Two grid cells, each with a fractional cloudiness of 25%. In the cell on the left, the 

cloud has the form of a rectangular block. In the cell on the right, the cloud is divided into four 
blocks. The cloud perimeter is twice as large in the cell on the left. It could be made even larger 
by dividing the same cloud area into a larger number of smaller cells. 

The perimeter of an individual cloud is 4 N/n, and the total perimeter of all clouds is 4 nN . 
Suppose that the cloud area increases through a process in which the areas of the individual 
clouds remain constant, while more clouds (of the same size as those already present) are added; 
in other words, some cells of the checkerboard are converted from clear to cloudy. Then if the 
cloud area is small, the total cloud perimeter increases in proportion to the number of clouds, 
which means that it increases in proportion to the total area occupied by the clouds: 

Perimeter ~ Area . (86) 

This scenario is analyzed in more detail below. 


5 Evaporation and/or sublimation can still occur even if a grid cell contains no clear air, e.g. due to large- 
scale sinking motion, but such phase changes do not necessarily reduce the cloud area until the final and 
complete destruction of the cloud occurs. 
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As a second example, suppose that the cloudy portion of the grid cell is arranged in a 

single square block of area N 2 , as illustrated in the right-hand panel of Fig. 3. In this case, the 
perimeter of the cloud is 4 N , so that the total cloud perimeter varies proportion to the square root 
of the area: 


perimeter ~ VArea . (87) 

We now analyze the first possibility above, in some detail. Again, consider a checkerboard 
cloud field, consisting of an array of square cells, each with area a . Each cell can be either cloudy 
or clear. Suppose that a certain number of cloudy cells are randomly arranged in the field, while 
the remaining cells are clear. If we randomly choose a clear cell to be converted into a cloudy cell, 
what happens to the total area and total perimeter of the cloud field? The total cloud area increases 
by a , but the change in perimeter depends on whether or not the new cloudy cell has cloudy 
neighbors. If each of the four neighbors (across cell walls) is clear, then the perimeter increases by 

Aja . If one of the neighbors is already cloudy, then the net change in perimeter is ija . If two 
neighbors are cloudy, there is no net change in the perimeter. If three are cloudy, the perimeter 

actually decreases; the net change in this case is -2 Ja. Finally, if all four neighbors are cloudy, 
the net change in the perimeter is -4 Ja. For a random arrangement of clouds with fractional area 
/, the probabilities of zero, one, two, three, and four cloudy neighbors are as indicated in Table 1 
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and Fig. 4. A simple calculation shows that the expected change in the total perimeter, when one 


Number of 
cloudy 
neighbors 

Probability 

Change in 
perimeter 

0 

/ 

4«/a 

1 

4/(1-/) 

2*Ja 

2 

«/(l-/) 2 

0 

3 

4/(1 -/) 3 

-ija 

4 

(l-/) 4 

-4ja 


Table 1 : Probability and incremental change in total cloud perimeter for an 
array of square cells. Here neighbors are counted only across cell walls. 
See Fig. 4. 
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Figure 4: Here we display all possible combinations of square cells with two neighbors across 

cell walls; the positions of the neighbors are indicated by the dots. Duplicates are shaded; there 
are 6 unique arrangements. This is where the factor of 6 comes from in the third row of Table 1. 


cloudy cell is added, is 
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d(perimeter) = 4 Ja( 1-2 /) . 


( 88 ) 


Integration gives 


perimeter(/) = 4jaf{\ -/). 


(89) 


The perimeter thus increases most rapidly with / when / is small, reaches a maximum when 
/ = 1 / 2 , and decreases symmetrically with further increases in /. 

Obviously the preceding example, with its array of square clouds, is not very realistic. 
Real cloud fields do sometimes resemble arrays of hexagonal cells, however. In contrast to square 
arrays, hexagonal arrays have the nice property that every neighbor lies directly across a cell wall. 
For a hexagonal cloud field, we obtain the probabilities and perimeter changes shown in Table 2 
and Fig. 5. The result is almost the same as before: 


Number of 
cloudy 
neighbors 

Probability 

Change in 
perimeter 

0 

/ 

61 

1 

6/(1 -/) 

41 

2 

15/ 4 (1-/) 2 

21 

3 

20/(1 -/) 3 

0 

4 

15 /d -/) 4 

- 21 

5 

6/(1 -f) 5 

-41 

6 

(l-/) 6 

-61 


Table 2: Probability and incremental change in total cloud perimeter for an 
array of hexagonal cells. Here l is the length of one segment of a 
hexagonal cell wall. See Fig. 5. 
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Tab!® 2. Tfa® factor of 2® in the fourth row earn be obtained b^ requiring that the sum of ail 
probabilities equal® 1. 


«/(perimeter) = 6/(1 -2f) , 


(90) 


perimeterf/) = 6//(l-/) ; 


(91) 


here / is the length of one cell wall. 
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4.2 Parameterization 

Assuming that the cloud field is divided into an array of cloudy cells, we expect the total 
perimeter of the cloud field to vary with the total area of the cloud field in accord with the 
/(!-/) dependence shown in (89) and (91). This motivates the following assumptions: 


•E^clcfclr C c\d.c\i m A 


cld,clr 


A cldA ck] 

(q 1 + q ■ \ 

evap evap 

1 

<N 

/ N 

i , i 

q c+q i ) 


(92) 


i A 


E clr,cld = C clr cld m A 


A c\dA clr 

. 1 , . 1 
Qc cond ^ Qi cond 

i 2 

(A 1 ) J 

Xq c + qi i)( 1 — RH* clr). 


(93) 


Here C cldclr and C clrcld are nondimensional parameters, assumed for simplicity and in the 
absence of evidence to the contrary to be constants; (<7 c )* evap ^ 0 and (# z ) evap ^ 0 are the rates 
of change of q c and q t due to evaporation in the cloudy part of grid cell i ; (9 c )* cond ^ 0 and 
(9))*cond - 0 are the corresponding rates of change of q c and q i due to condensation in the 

cloudy part of grid cell i ; and RH‘ c h = {q v)ck/ {q * )ck is the clear-air relative humidity of grid 

i i i 1 

cell i . Both E cld,clr and E dr.cld are assumed to be proportional to A c ld A c lr , which, in turn, is 

assumed to be proportional to the total perimeter of the cloud, based on the arguments given 

above. According to (92), evaporation is associated with a transformation of cloudy air into clear 

air, which becomes very efficient when q c + q i -> 0 . Similarly, (93) states that condensation is 
associated with a transformation of clear air into cloudy air, which becomes very efficient when 

RH 1 clr — > 1 , i.e. when the relative humidity of the clear air approaches 100%. When the relative 
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humidity of the clear air reaches 100%, the denominator of (93) goes to zero, which means that 

E l c ]i,c\d can become non-zero even if A l c id = 0; this means that we can create cloud directly 
from clear air as the water vapor mixing ratio of the clear air reaches saturation. 

We can develop a further interpretation of (92), as follows. Begin by writing 


— [m 1 ^ cld (qc + qi)] 


1 evap 


= m 


^cld[(<7 c )‘ evap + W.O'evap] + ^ c + * A dd) ] evap 

= /M , ^dd[(^ c )' evap + (9,) I e vap] - ^ Ic + ^ i ')- £:cld ’ clr • 


(94) 


The left-hand side of (94) represents the total rate of evaporation of cloud mass inside the grid 
cell. According to (94), this total consists of a part that comes from the reduction of the cloud 

water mixing ratio within the cloudy region, given by m‘ A c \d[(q c )‘ eySLp + (^O’evapl > 311(1 3 P 311 

i i i 

that comes from a reduction of the cloud area, given by -{q c + q i)E cld,clr • Let the ratio of these 
two contributions be denoted by i? evap , i.e. 

_ rate of cloud destruction by transformation of cloudy air to clear air 
evap " rate of cloud destruction by reduction of the in-cloud condensate mixing ratio 

i i « (95) 

_ (<7 c + i)E cld,clr 

;4cld[(<7 c ) evap (?/) evap^ 


Rearranging (95) gives 

q 1 +q i \ 

“c evap “i evap 

v q c + q i J 


E l dd,clr = -5 evap mU c ld 


(96) 
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Comparing (96) with (92), we see that 


_ rate of cloud destruction by transformation of cloudy air to clear air 
eva P rate of cloud destruction by reduction of the in-cloud condensate mixing ratio 

.i 

r — clr 
" C cld,clr f • 

A 


(97) 


Recalling our assumption that C cld clr is a constant, we see that (97) states that when A 1 c ir/ A 1 is 

small, i.e. when there is very little clear air in the grid cell, the destruction of the cloud is mostly 
due to a reduction in the in-cloud condensate mixing ratio, and reduction in the cloud area is 
secondary. In effect, we are assuming that the transformation of cloudy air into clear air is favored 
by pre-existing clear air. This assumption implies that, if the grid cell is completely filed with 
cloud, clear air can be produced only by reducing the condensate mixing ratio to zero, at which 

point £ ! c id,clr -> oo . As the area covered by clear air increases, the transformation of cloudy air to 
clear air becomes increasingly important, relative to the reduction of the in-cloud condensate 
mixing ratio. 

To similarly interpret (93), begin with 

iim'A cld(9c + ?'.')]} 

dl J cond 

- + (?,.) i oond ] + (?'c + « i ')[|<'"^'" d )] cond (98> 

- m , ^cld[(? ( .) l codd + (? / ) l cond ]-(? , c + 9'i)£ clr,cld , 


and define 
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B 


_ cloud generation by transformation of clear air to cloudy air 
cond = c ioud generation by an increase of the in-cloud condensate mixing ratio 

_ (<7 c + q i)E clr,cld 


( 99 ) 


m ^cld[(<? c ) cond + cond-1 


so that 


E clr,cld ~ ^cond m ^ cld 


' (g c V cond + (Aji cond 
(?'c + 9i) 


( 100 ) 


Comparing (100) and (93), we see that 


B 


cond 


cloud generation by transformation of clear air to cloudy air 

cloud generation by an increase of the in-cloud condensate mixing ratio 


= C 


A clr 

clr,cld i 

A 


1 — RH 1 clr 


(101) 


An interpretation of (101) is that, for a given value of A'ck/A 1 , as RH l c lr — > 1 cloud generation 
by transformation of clear air to cloudy air dominates over cloud generation by an increase of the 
in-cloud condensate mixing ratio. In other words, transformation of clear air to cloudy air is 
favored if there is lots of humid clear air available. On the other hand, as the cloud amount 

approaches 100%, for a given value of RH 1 c ir, cloud generation by an increase of the in-cloud 
mixing ratio dominates over cloud generation by transformation of clear air to cloudy air. 

It is useful to compare (92) and (93) with the corresponding parameterizations of Tiedtke 
(1993), which are (using the notation of the present paper) 


E cld,clr — 0 , 


( 102 ) 
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fi'cMd - — hr] M ‘'4‘S*’° 1 • 

j q'H-RrfJ V dt J 


( 103 ) 


Here AH * is (in Tiedtke’s parameterization) proportional to the rate of production of cloud water 
dt 

and/or cloud ice by condensation. In contrast to Tiedtke (1993), we permit a non-zero 
“evaporative mass flux,” £ ! c ld,clr • On the other hand, our parameterization of the “condensation 

mass flux,” £' c ircld is similar to Tiedtke’s in that both are proportional to - — — 1 

A‘ M — Rn clr 

times a measure of the condensation rate. 

5. Computational procedure 

Using (81)-(85) and (92)-(93), we now rewrite (73)-(80) as follows: 


l 


O , l A x _ ~ l > 1 l ’> 1 i 

— (/w A clr) “ ~/ M ^ nl 


U t U f ?A ’ clr 


1 A 1 

-m A 


A cldA clr 
i 2 

L (A 1 ) J 


'cld,clr 


-A, 1 

A 


f q 1 + q 1 ^ 

evap ^ i evap 

k qc+qi J 


+ C 


clr,cld 


. i + . / 
cond _ £/ cond 

hqc + qiW-Rrfch). 


(104) 


— clr 3/* — A* M meso) + 

oz 


, , id A clr ~.iA clr^ 

M„ -s— —D — 

c dz 4 i 
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d , iJ s x— > * /, j, ? J,fA cld 

^4 dd) = - 2 ,™ v «* 


+ m l A l 


A cld A clr 

ic r 

i 2 

(A‘) J 

j L clr, cld 


?c cond ~ ^ i cond 


A 

+ c 


l, l 


( , 


. I 


< c evap + evap 
cld,clr| i i 

\ q c+q i ) 


-^-(A^ldM* + A 1 M meso) + 
dz 


, , idA l c\A. 

M -r— + 


D* 

c dz 

l A 1 J 



O , i .i. a i, C i, i' ji, f d , T7 1 A\ 

—(m A ) = v nl -j-(M A ), 


- i, r . , 

d, i,i J \ _ ,i,f i,? jiJA clr ?> 1 , 

-z~( m h clr-'d clr) — ~z_J n v n ^ T~F~ ’ l 

Ot , y>‘ 

i A 


+ m l A l 


A q\&A clr 


fq ' +q i \ 

“ c evap “ / evap 

r 

4 c cond + 4\ cond ] 

- U) . 

Q 

c cld, clr 

i , i 

V q c + q i J 

^cl^cld 

.(q‘c + q i i){ 1 - RH* clr) J J 


-^hA^lrM meso)// clr ] + (S^^dr 


+ 


ckA* clr) ~ D‘h‘clx^-f 

dz A 1 J 


* i, i . , 

O , i,i J s v- 1 * »> *. ^ *’-'4 cld ' ,, 

-=-(/n h ddA dd) = ->w v «/ — 77 -fl cld 

a/ r y 




A q\&A clr 

| ^clr,cld 

r • i + • i 

4 C cond 4j cond 

c 

(a 1 + q! M 

evap evap 1 

- (A) 2 - 

Xqc + qiW-RH* clr). 

C cld,clr 

«• i i 

^ q c + q i Ji 


-•^(A'cldM + A l M meso)h cld j + (5 , A ) ! cld ^ J c ld 


+ 


M c ’ cld A dd) d" D 


.A 


if , i ,i A cld 
h cu — h cld — 7 


V 


J J 


(105) 


(106) 


(107) 


( 108 ) 
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( 109 ) 


!-( nth 1 A ') = -Y m’ l v ,l nl l,l h ,, ‘ + Yc'’'(A) 

dt i i 

-^{^4 [A4*h + A/^mesc^A eld - A clr)]} + M cld(<S , ^) c ^ + ^4 clr(*S'/,) c | r ] 

+ A 1 Me $■ + Dtyca ~ h l ) + e\ — (A* - /t'clr) + h l - Add) 

5z ^ ^ 


.M" 


Following Tiedtke (1993), we prognostically determine ^ f c ld . In view of (66), there is no 
need to predict A l c lr and ^4* c ld separately; prediction of ^4 z c ld suffices, since (66) can then be 
used to diagnose a‘ c \t- 

Consider two possible ways to use equations (107)-(109): 

Method A: Separately predict the properties of the clear and cloudy sub-regions. To 

determine Ji , we could use (67). This approach is simple, and it is fine in 
principle, but practical and/or philosophical objections can be raised, along the 

following lines: Existing models predict V . Rewriting a model to predict A* c lr 

and A*cld would be a lot of work. Other current and/or future cloud 
parameterizations do/will not make use of the framework proposed in the 
present paper, so building this framework into a model at the “foundation level” 
seems unwise. 

Method B: Continue to predict A*, as in current models, using (109). In addition, predict 

A l c id and A^id - h‘ck, the latter using an equation to be derived below. From 
these predicted values, diagnose A ! c id and A*clr using 
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( 110 ) 


h' c Id = h' + 


1 - 


A cld^ 


( h cld — h clr) , 


and 


1 1 yi ^ clfl / 7 i t ^ \ 

h c\i ~ h — (^ cld — h clr) 5 

A 1 


(111) 


respectively. 

We use Method B to predict the generalized moist static energy and the water vapor mixing ratio. 
We use Method A to predict the condensed water variables. Further discussion is given in Section 
6 . 


For use with Method B, we need a prognostic equation for /i'cid-^clr, which can be 
derived as follows: By using (104) in (107), we can obtain the advective form corresponding to 
(107): 


dh clr _ 1 y 


h f . \ 


f, r U f / U ? A clr 


* i, i hi 1 
m v n l 


A 1 ’ 1 ) 


(h clr-Zl'clr) 


-C 


^cldV 


cld,clr' 


V ^ j 


^ \ 


j c evap ^ 7 ; evap 

V + <7*/ 7 


(^*cld - /^clr) 



— (A/cld-M'clr)!] 

1 ^ 

Im m 

L A 1 JJ 


[ 5/z clr clr ^_c_dh clr 

| ~dz _j „* <3z 




m 


( 112 ) 


Here we have used (45). Similarly, we can show that 
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dh cld 1 


_ cia — 1 y 

dt i A , 

m ^4 cld i 


f \ h * \ 

a /, f /, r cid 

m v n l 




~j, * 


, f h t ji x 

( h cld — h cld) 


J 


-C 




clr,cld 


V A J 


. 2 , . I 

cond cond 


(A cld “ h clr) 


-K + i 


'^clr 


2 2 


(Af* cld — A/* clr) 


m‘ ml ^ 


L(^c + ^/)(l-^clr) 

5A Z cld , cld , M c dh 1 cld , £>' 


( 113 ) 


fdz 


m 


m 


J dz 


+ 


m A l c Id 


(h cu-h cld) • 


In the horizontal advection term of (1 12) and in the horizontal advection and cumulus detrainment 

terms of (113), there is an apparent danger of division by zero when ^ c lr -> 0 or A l c \ d -> 0 . The 

physical meaning is that h l c ]x can be instantly “re-set” by advection when ^ c ir = 0, and that 

A 1 cld can be instantly re-set by either advection or cumulus detrainment when A 1 c ld = 0 . The 
danger of division by zero can be eliminated through the use of suitable finite-difference methods. 
Details are discussed in Appendix B. 

To obtain a prognostic equation for h‘c\d -n c h, subtract (112) from (113): 
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-~-{h c\A — h c\i) ^ 


dt 


( *i,f i,f A?*/| 

mvrtl__A cld 

m A 1 c Id 


(V’* cld -Veld) 


+ i 
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( A U i U ? jU i y> f t ^ 
m v w / ^ clr 

m^clr 


z ? i, i ji N 

(/* clr — h clr) 


^ eld ~ ^ clr ^ 

i , i , 
V q c + q i ) 


( A clrlf 9c cond + 9/ cond 


'clr,cld 




1 — RH c\x 


-C 


eld, clr 


V C kT 

V iV > 


(<T + <V 

V ^c evap ^ i evap 


) (H4) 


M* d , , i ,i v 

:■=-(« cld — h clr)- 

( (7Z 
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(■ _ M c\x ¥ v4*cir dVcld + Vcld dVdV 

, ' » <9z <9z 

V m m J\ A A ) 


,(Vcld Weir, 


m 


/w 


m; a , 

c o , , i 7 1 


m 


‘ dz 


(h cld- A clr) + 


^ D ; ^ 

V/V.Vcldy' 


(A cu — h cld) 


The lateral mass exchange terms of (1 14) act to damp (Veld - Vclr) towards zero. 

In order to predict .Veld, V and (Veld - Vclr), we must determine aV meso 9 Mcid',^ 
M clr * ■ Methods to do so are discussed in Section 7. 


6. Microphysical processes 

The microphysical parameterizations of the model, which represent the microphysical 
processes occurring in the cloudy portion of the grid cell, basically follow those described by 
Fowler et al. (1996). We distinguish a total of five prognostic water species: water vapor, with 
mixing ratio q v ; cloud water, with mixing ratio q c ; cloud ice, with mixing ratio q { \ rain falling 

from stratiform clouds, with mixing ratio q r ; and snow falling from stratiform clouds, with 

mixing ratio q s . Our bulk cloud microphysics parameterization also includes a prognostic 
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equation for the generalized moist static energy, h, defined by Lord (1978) and discussed by 
Fowler et al. (1996). In addition, of course, we prognose A* c l d . 


Each of the six thermodynamic variables (h, q v , q c , q ( , q r , and q s ) is assigned values 

for the cloudy and clear regions separately, although as discussed in Section 2 we define the clear 
sub-region to be one in which the mixing ratios of cloud water and cloud ice are zero. Table 3 
summarizes the prognostic thermodynamic variables used by Eauliq NG. The model prognoses 


Variable Definition 

Symbol 

Mean generalized moist static energy 

h 

Generalized moist static energy difference 
between cloudy and clear air 

h c\d ~ h cb 

Mean water vapor mixing ratio 

% 

Water vapor mixing ratio difference 
between cloudy and clear air 

^v) c id - (^v) c lr 

Stratiform cloud area 

^cld 

Cloud water mixing ratio in the cloudy 
sub-region 

Vc’ 

The subscript “cld” is 
omitted on the grounds that 
it would be redundant. 

Cloud ice mixing ratio in the cloudy sub- 
region 

<7,-; 

The subscript “cld” is 
omitted on die grounds that 
it would be redundant. 

Rain water mixing ratio in the cloudy sub- 
region 

^r) c ld 

Rain water mixing ratio in the clear sub- 
region 

<«r) * 


Table 3: Prognostic thermodynamic variables used by Eauliq NG. 
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Variable Definition 

Symbol 

Snow mixing ratio in the cloudy sub- 
region 

^) c ld 

Snow mixing ratio in the clear sub-region 



Table 3: Prognostic thermodynamic variables used by Eauliq NG. 


(i.e. time-steps) the grid-cell averaged generalized moist static energy, and the difference in 
generalized moist static energy between the cloudy and clear portions of the box. Similarly, the 
model prognoses the grid-cell averaged water vapor mixing ratio, and the difference in water 
vapor mixing ratio between the cloudy and clear portions of the box. In the terminology of 
Section 5, we use Method B for the generalized moist static energy and the water vapor. 

The model also prognoses the cloud water mixing ratio in the cloudy portion of the box 
(which can be considered to be the difference between the cloud water mixing ratio in the cloudy 
and clear portions of the cell). Similarly, the model prognoses the cloud ice mixing ratio in the 
cloudy portion of the box. There is no need to separately prognose the grid-cell averaged values of 
the cloud water or cloud ice, but we can diagnose them, if we wish, using (67). This means that 
Method A is used for the cloud water and cloud ice. 

Precipitation obviously originates in the cloudy portions of grid cells, but it can fall into 
either the cloudy portion or the clear portion of a lower-level cell. We use Method A to predict the 
rain water and snow mixing ratios in the cloudy and clear sub-regions of each cell. 

7. Mesoscale circulations 

7. 1 Background 

As discussed in the Introduction, we assume that the stratiform cloud under consideration 
is and remains through time neutrally buoyant with respect to its clear environment. Neglecting 
(just temporarily) virtual temperature effects, this neutral buoyancy condition can be expressed 
by: 
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(115) 


T = 


T c\i T c\i : 


r 


“ r cld ” r clr’ 


(116) 


3T 

where r = — r- is the lapse rate of temperature; and finally 
oz 


d_j, _ ^cld _ ^clr 
dt dt dt 


(117) 


We cannot simply assume that (115)-(117) are satisfied; a physical process must act to ensure that 
they are satisfied, and such a process can produce additional effects that go beyond (115)-(117). 
A candidate process is proposed below. 

Consider a non-precipitating stratiform liquid water cloud that is caught up in dynamically 
imposed large-scale vertical motion, which could be either upward or downward. For simplicity, 
suppose that radiative cooling and horizontal advection are negligible, and neglect cumulus 
effects. We temporarily assume that within the cloud layer the cloud amount is 100%. As the 
cloudy air rises or sinks, it cools or warms along a moist adiabat, so that 


DT _ 

— — = -wr in cloud, 
Dt m 


(118) 


where the Lagrangian derivative is defined to follow the mean motion: 


aj.|( ) + V.V( ) + 4< ). 


(119) 
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In (119), which is (for purposes of this section only) the definition of the Lagrangian derivative, 
we include a term that represents horizontal advection, even though we are currently assuming 
that horizontal advection is negligible. We have used the mean vertical velocity, w , in (1 19), with 
the understanding that because the cloud amount is 100% by assumption, w = w dd . 

The dry air above cloud top follows a dry adiabat, so that 

DJi - -wT, in clear air. (120) 

Dt d 

Here we have used the mean vertical velocity, with the understanding that because we are 
considering a level above the cloud top, vv = w clr . Suppose that at a given moment the 

temperature varies continuously with height. For w < 0 , the sinking air will subsequently warm 
both above and below cloud top, but more rapidly (along the dry adiabat) above cloud top, and 
more slowly (along the moist adiabat) below cloud top, so that after some time there will be a 
discontinuous upward increase of temperature across the cloud top. This illustrates that large- 
scale sinking motion tends to produce sharp cloud-top inversions, which are of course often seen 
at real cloud tops. Although there are additional mechanisms, such as radiative cooling, that can 
promote the formation of inversions at stratiform cloud tops, we note here that large-scale sinking 
motion alone suffices. 

For w > 0 , the rising air cools both above and below cloud top but more rapidly above, so 
that an initially continuous temperature profile tends to develop a statically unstable upward 
decrease of temperature across the cloud top. Of course, small-scale convection will very quickly 
develop so as to prevent such a statically unstable layer from actually being generated (Arakawa 
and Schubert, 1974). 
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cld 

Now allow the possibility of fractional cloudiness, denoted by / = — — . We continue to 

A 


suppose for simplicity that radiative cooling, horizontal advection, and cumulus convection are 
negligible. Assume that the temperature distribution is initially horizontally uniform and 
vertically continuous. Consider the case of large-scale subsidence. As the sinking cloud 
evaporates, its temperature increases along a moist adiabat. If the environment at the same level as 
the cloud sinks at the same rate as the cloud, its temperature will increase more rapidly, along the 
dry adiabat. This means that after some time the cloud will be colder than its environment at the 
same level, so that (11 5)-(l 17) will be violated. We can avoid this predicament if we allow the 
vertical velocities experienced by the cloud and the environment to be different, i.e. 


dT 


cld 


dt 


w cld( r cld -r m)’ 


(121) 


gJ Clr ■ M 'clr( r clr- r d)- <> 22 ) 

The large-scale vertical velocity is given by the area average of the vertical velocities of the cloud 
and its environment: 


w = /^dd + O-ZKlr- 


(123) 


Equations very similar to (121)-(123) were studied by Bjerknes (1938), although he did not use 
(1 15)-(1 17). Combining (1 15)-(1 17) and (121)-(123), we find that 


w 


clr " 


w = - 


/(r rf -r M ) 

/(Trf-n-O -f)(T-r m ) 


w 


(124) 


and 
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w 


(125) 


- = O-ZK^-rj 
Wcld w [/(r^-r)-(i -/)(r-r m ) 

are the mesoscale vertical velocities needed to keep the cloud and environment at the same 
temperature. From (124) and (125), we see that = 0 if w = 0; this is natural, 

because in the present simplified scenario it is only the large-scale vertical motion that is 
disturbing the equilibrium of the system. Suppose that T = T m , which will be the case if 
turbulence keeps the interior of the cloud in a well-mixed state. Then ( 1 24) and ( 1 25) reduce to 

w c lr = 0 ’ (126) 


w 

W dd " /' 


(127) 


These results mean that, within the current but temporary restrictive assumptions, vertical motion 
occurs entirely within the cloudy region, and not at all in the environment. The cloudy air rises, or 
the cloudy air sinks. The environment sits still. If a large-scale dynamical process imposes large- 
scale vertical motion on a region containing fractional cloudiness, the cloudy air carries out the 
required vertical movements, depending on the value of /, but the environment does not join in. 
Such differences in vertical motion between the cloud and its environment can be interpreted as 
consequences of an adjustment process, which counteracts the warming or cooling of the cloud 
relative to its environment so as to maintain a state of balance, with no horizontal pressure 
differences between the stratiform cloud and its clear environment. 

As discussed earlier, the mesoscale circulations represented by (124)-(125) or (126)-(127) 
transport various quantities vertically. From (124) and (125), the mesoscale mass flux can be 
written as 
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(128) 


w cld- w clr = 


/(r rf -r)-(i-/)(r-r m ) 


w . 


In case T = r , this reduces to 
m 7 


w 

w cld — w clr = 


(129) 


which implies that, for a given w, a strong mesoscale mass flux is favored by a small cloud 
fractional area. If the cloud amount approaches one, the mesoscale mass flux goes to zero. Note 
that the mesoscale mass flux has the same sign as the large-scale vertical motion; in the case of 
large-scale subsidence, the mesoscale mass flux is downward. 

7.2 Parameterization of the mesoscale mass flux 

The preceding discussion is intended to show that the mesoscale circulations should act in 
such a way as to remove buoyancy differences between the stratiform cloudy and clear air. At 
present, we parameterize the mesoscale mass flux by introducing a simple prognostic equation for 

w dd - u/ c lr, based on (1 14), which allows w' c ld - wclr to increase due to mesoscale buoyancy 
forces: 
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cld - w'clr) = 
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( *ij ij f -J, ? \ 
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1 uz 
m 


clr J 


(Mj M^a 


\ 

^ddl 

C cld, clr 

J 

a' j 


evap + evap) \ 


i i 

V m m J 


A clr^vv cld cld^w clr 


.J dz 


A dz 


-|[KS] cld -[Kp)]Ji [<r ^- (r ^ 1 


(130) 


+ 


Ka, i 


m 


i dz 


(w cld — w clr) + 


D 


i ^ 


w A l c \&) 


(w'cu-w'dd) 


Here T v is the virtual temperature. We include condensate loading in the definition of T v . The 
pressure terms of ( 1 30) must be parameterized or neglected; we currently neglect them. We expect 

that the effect of w'dd-w^lr, as predicted using (130), will be to limit (T v ) - (^ v ) c lr t0 

relatively small values — relative, that is, to what might occur if w^ld - w* c ir = 0 . In other words, 

we expect that the use of (130) will produce values of w^id-wW which act to keep the 
stratiform cloud neutrally buoyant with respect to its clear environment. 

Once w ! c id - vv^lr has been determined using (130), we can obtain A4 * meso from 


M^cld-A^clr = m l (w l c ld-w c ]r) 


(131) 
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with (45), and and M c ^‘ from 


M dd = m + 1 “ 


^cld* 


(M* cld — M* clr) > 


( 132 ) 


and 


A/ clr = V-— (A/dd-A/'clr), 

A 1 


( 133 ) 


respectively. We assume that M* is determined by the large-scale model. 

8. Summary and conclusions 

We have formulated a generalized version of Eauliq, the stratiform cloud parameterization 
developed by Fowler et al. (1996). The new version, called Eauliq NG, has the following 
enhancements: 

• A prognostic stratiform cloudiness, similar to that of Tiedtke ( 1 993). 

• Separate prognosis of the thermodynamic properties of the cloudy and clear air. 

• Fully consistent incorporation of the cumulus terms as they affect the microphysical 
variables and the cloud amount. 

• Diagnosis of distinct vertical motion fields for the cloudy and clear portions of each 
grid cell, formulated in such a way that the stratiform cloud tends to remain neutrally 
buoyant through time. These distinct vertical motion fields are used consistently in the 
prediction of all thermodynamic fields. 
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We have not considered mesoscale fluctuations of the horizontal velocity, even though 
they are implied (through continuity) by the mesoscale fluctuations of the vertical velocity which 
we have included in our model. This deficiency should be remedied in the future. 

This technical report is being published in order to record the preliminary formulation of 
Eauliq NG prior to its testing in a single column model and subsequently in a general circulation 
model. The final version of Eauliq NG may depart somewhat from the description given here, 
depending on the outcome of such tests. 
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Appendix A 

Analysis of Vertical Fluxes 


We can rewrite (46) as follows: 


A ckAd ck ~ A Ad* — A c\dAd cld ~ A cu M* CU 5 


or 


A M clr — A M* — A cldiM* cld ~ clr) A cu(A^cu — A/*clr ) * 


or 


i j 


JckMck = A* c\rM - A ° ldA Clf (M J cld - Af'clr) - ^ (M^u - M*ck) 


A clr A 


Similarly, we can rewrite (47) as 


A c\dAd cld — A Ad* — A clr Ad clr ~ A cu Ad 1 cu , 


or 


A Ad cld — A Ad* + A clr (Ad* cld — Ad * clr) — A cu (Ad* cu — Af*cld) > 


or 


J J 


i 


A 1 cXdAd cld = A l c \dAd + - ■-^-■- - (i^cld-i^clr) ^ T ^(Ad i C u- Ad 1 M ) . 

A A 


( 134 ) 


(135) 


(136) 


(137) 


(138) 


( 139 ) 
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The large-scale flux of h is given by 


Al mwh = c\xM clr^clr + -^cld-W cld^ cld + A cu^cu^ cu 


Substituting from (136) and (139), we find that 


(140) 
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A l mwh = 


l A 


i t~J AcldAckstJ 


A'ckM* ~ 


(M * dd — M 1 clr) — A clr — cu — A/ 1 clr) 


h clr 


( 141 ) 


i A 


A'cldM* + ^4 clr 


(A/dd - A/ clr) - ^ dd— (A^cu - M* old) 


h l c\d + Al CU M 1 cu h 


cu 


= M* (A 1 c \rh l c \x + ^cld^'cld) + ~ ^4 — (M^dd - A/clrX^dd - h l c \ r) 


- ^'clr— (A^cu - M* dr^'clr - ^cld— (A /cu - A/dd)A'cld + A'—a/cu^cu 

A 1 A 1 A 1 

= M* (A‘ dr h‘c\i + A'ddh'cld + cu) + - — - A^ c lr)(Add - A dr) 


H r[ — -4 clrCA^cu — M* clr) A clr — A cld (A/* cu — A/^ cld) A cld + -4 (A/^ cu — M*)h cu] 
A 1 

= AlMIi + - cld f Clr (A/cld - A/clr)(A , dd - A'dr) 


+ — 7" { A/^cut - ^ clr^ clr — A c\dh cld + A h cu ] 

A 1 


+ A ofoM 1 ckh clr + A cld^^cld h cld — A Af* h cu} 


= A l MV + - ° ldA — (M* cld ~ A/clr)(/i'cld - A'clr) 


+ Afcul-A'h 1 + A'cuh'cu + ^Vcu] + (A'mV - / cuM* C uh l C u) - aWIi'cu 


— A M h d — — (A^cld — A/^clrXA cld — h clr) + A cuCAf* cu — A/)(/t cu — h ) 

A 1 


= A\]^h l + A/mesoC/ldd - h‘ dr) + cu - A )] . 
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Appendix B 

Time-Differencing 


The prognostic equation for the cloud amount, (105), is repeated here for convenience, 
with a slight expansion and rearrangement of the terms: 

/) = Vcld + mVXl ~/)X i cld + (l C 142 ) 

St A l 

where for convenience we define 


a _ a cid 

/= A 1 ’ 


(143) 


Xcld = 


1 

n 

4c cond + cond 

(?c + qt) 

C clr,cld 

n 1 — RH* clr J 


^cld,clr^ 


i 

c evap 


+ ^ i evap^ 


(144) 


and 


id/ 

6z 

Note that ^cld can diverge to either plus or minus infinity. Here we discuss the time differencing 
of the lateral mass exchange and cumulus detrainment terms, which are shown explicitly in (142). 
We write 


A(/V + A/'meso) • 
oz 


(145) 


i 1 

V)/ cld = -— 2, 
A { 


* i, i 

lr.i. i 1 U t ,i, i'A cld ^ X/f 
m v n l — r- + M „ 


A 1 ’ 1 
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(146) 


(« , /)" +1 -(m'/)" _ / 

— V cld 

At 

*m ' ■(/)"' * ' [ 1 - (/)’ ' " ‘ ] X cld + [ 1 - (/)" * ' ]^. • 

A 


Here superscripts n and n + 1 denote successive time levels, and At is the time step. Eq. (146) 

n+ 1 


can be solved as a quadratic equation for (/*) 


[(/)” + 1 ] 2 »<VcldA/ + (/) 


n+ 1 




I 2 


Z> 


/ A 


1 - m x cldAr + — .A t 
A 1 ) 


(m /) + VcldA t + —.At 


D 


= 0 


(147) 


The solution is 


(/) 


n+ 1 


l-mxcld&t + —At 


2/«YcldAt 


f i . A 2 

1 - /n x cldAt + — -.At 

v Y J 


+ 4(m'x , cldA0 


2m i x , cldAt 


{m f) + ij/cldAt + — At 


(148) 


Consider a simplified case in which lateral mass exchange is the only process. Then (148) 
reduces to 


+ 1 _ 1 - mV cldA? 


(/) 


a/(1 ~/»VcldA t) +4(mVcldA;)(m*/) 


2m , x , cldAt 


2m'x ; cldAt 


(149) 
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For the special case 


— > oo this further simplifies to 


IX cld| 


j n + 1 1 1 

( f) = - ±- 
U ) 2 - 2 


j n + 1 J n + 

If we choose the plus sign, we get (/) = 1 ; if we choose the minus, we get (/) 

It thus appears that we should choose the plus sign for x*cld > 0 , and the minus for % c ld < 
Now consider a second special case in which %cld = 0 . Then (147) reduces to 


(/) 


i A n i D 

( m j) + v|/ cldAt H :A t 

«+l = A 1 

1 + —At 


n+ 1 


For D l ->• oo we get (f 1 ) = 1 , as expected. 


Next, consider (1 14), the prognostic equation for h l c id - ^clr : 


^(^'cld - h l c \x) = f(h)~ [i\h l dd- h l c \i) + 

ot 


D 


i \ 


\jnA l c \&) 


( h cu — h cld) • 


Here for convenience we define 


( 150 ) 

1 

= 1 . 
0. 


(151) 


(152) 
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( 153 ) 


A*)— Z 


/ A i, r i, r iU? ^ 

' m_v^jiV_A cld ' 

m^cld 


/i'’ 1 ,J \ 

(A cld - h dd) 


+ Z 

f 


7 - J, f i, ? ,i, i i 1 \ 

m v n l A clr 


rn A l c lr 


( h clr — h clr) t^-(/i cld - h clr) 

J ™ ° Z 


^cld ^clr 


. siJ 


V m m A A 




A clr 5 A cld A cld^/i clr 


( i dz 


,i dz 


+ 


Weld Wclr 


m 


m 


M 


m 


i dz 


(h l cld - /Air) 


H = 


(A + q l i) 


"clr, cld 


'A* v ' ' ■ • ‘ 

v ^ yv 


cond + co nd 

I - RH“a, 


\ 

( A 1 \ 

A cld 

c cld, clr 

) 

A 

\ A ) 


Wc'evap + i, ev,p> 1 < 1! *> 


Note that 


jx* >0. 


(155) 


We can rewrite (152) as 


dt 


{h cld — h clr) — 


A»+ 


D 


i \ 


\m A 1 c \&) 


(/l cu — h clr) 


1 f 


A + 


D 


i \ 


J V 




(Vcld-Vdr). (156) 


We write 


« + 1 


(A cld — h clr) — ( h cld — h clr) _ 


At 


A*)+ 


\.rn A* cldJ 


(h cu — h clr) 




/• D 

H + 


/ A 


TnA l cidJ 


(h cld — h clr) 


n + 1 


(157) 
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The solution is 


( h cld — h clr) + At 


i i n + 1 
( h cld - h clr) 


/(*)+ 


r Sl ^ 

VmA 1 c \d) 


( h cvi — h clr) 


1 +A t 


( i D ‘ ^ 

V + — 


i ji 


V 


m A cld ) 


( 158 ) 


For D l -> oo , we get 


(h cld — h clr) 


n+ 1 


= h 


cu ■ 


i'clr. 


( 159 ) 


For \x -> oo , we get 


( h cld — h clr) 


71+1 


= 0 . 


( 160 ) 
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